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Background: A signature of many dynamical models of dark energy is that they admit variation in the fine 
structure constant a over cosmological time scales. 

Purpose: We reconsider the analysis of the sensitivity of neutron resonance energies Ei to changes in a with a 
view to resolving uncertainties that plague earlier treatments. 

Methods: We point out that, with more appropriate choices of nuclear parameters, the standard estimate (due 
to Damour and Dyson) of the sensitivity for resonances in ^®°Sm is increased by a factor of 2.5. We go on to 
identify and compute excitation, Coulomb and deformation corrections. To this end, we use deformed Fermi 
density distributions fitted to the output of Hartee-Fock (HF) + BCS calculations (with both the SLy4 and 
SkM* Skyrme functionals), the energetics of the surface diffuseness of nuclei, and thermal properties of their 
deformation. We also invoke the eigenstate thermalization hypothesis, performing the requisite microcanonical 
averages with two phenomenological level densities which, via the leptodermous expansion of the level density 
parameter, include the effect of increased surface diffuseness. Theoretical uncertainties are assessed with the 
inter-model prescription of Dobaczewski et al. [J. Phys. G: Nucl. Part. Phys. 41, 074001 (2014)]. 

Results: The corrections diminish the revised ^®°Sm sensitivity but not by more than 25%. Subject to a weak 
and testable restriction on the change in m^/A (relative to the change in a) since the time when the Oklo 
reactors were active (m^ is the average of the u and d current quark masses, and A is the mass scale of quantum 
chromodynamics), we deduce that |aokio — Q!now| < 1.1 x 10~®Q:now (95% confidence level). The corresponding 
bound on the present-day time variation of a is tighter than the best limit to date from atomic clock experiments. 

Conclusions: The order of magnitude of our Oklo bound on changes in a is reliable. It is one order of magnitude 
lower than the Oklo-based bound most commonly adopted in earlier attempts to identify phenomenologically 
successful models of a-variation. 

PACS numbers: 06.20.Jr, 21.10.Ft, 24.10.Pa, 24.60.Lz 

Quite apart from their intrinsic metrological inter¬ 
est, studies of the possible variation of fundamental di¬ 
mensionless parameters like the fine structure constant 
probe beyond the standard models of elementary par¬ 
ticle physics and cosmology 0 , 0 . The initial applica¬ 
tion of the many-multiplet method to 128 Keck/HIRES 
quasar absorption systems suggested that, in the red- 
shift range 0.2 < z < 3.7, a was smaller than today by 
about 6 parts per million (ppm) 0. Subsequent analysis 
of Keck/HIRES and VLT/UVES spectra has lead to a 
refinement of this earlier claim, namely, that a appears 
to vary spatially across the sky 0. An angular dipole 
model of amplitude ^ 10 ppm is favoured at the d.lcr 
level over a simple monopole model in which a does not 
change across the sky but could be different from cur¬ 
rent laboratory values; the data also supports a dipole 
model with an amplitude proportional to the look-back 
time. In view of the paradigm-shifting ramifications of 
these hndings, mined from data taken for other purposes, 
and concerns about the wavelength calibration of the 
HIRES and UVES spectrographs (which seem to have 
been borne out by the recent identification of long range 
distortions i), the UVES LARGE Program of dedicated 
observations has been initiated to check the evidence for 
non-zero changes in a 0 . However, the difficulties of the 
spectroscopic measurements involved demand a new gen¬ 
eration of ultra-stable high resolution spectrographs 0 , 
two of which (PEPSI and ESPRESSO) it is envisaged 


will begin operating within a year or so. 

The reconciliation of these claimed changes in a with 
stringent bounds from the Oklo natural fission reac- 
and single ion optical clocks has pro- 


ill 


tors 

voked considerable theoretical effort (work prior to 2011 
is reviewed in Ref. 0 ), with several models under ac¬ 
tive development |l3l4l6| . Nevertheless, there is some 
ambivalence in the literature about the importance of 
the Oklo geochemical data. Significantly, the Oklo phe- 
nomenon is completely ignored in recent studies |l7| of 
the potential of high precision measurements of the red- 
shift dependence of a (among other quantities) to dis¬ 
tinguish between different dynamical dark energy mod¬ 
els and pin down unihcation scenarios. Other recent 
papers dismiss Oklo bounds for being “strongly model 
dependent and possibly subject to criticism” [l6|, “sub¬ 


ject to much larger theoretical and systematic uncertain¬ 
ties than ... spectroscopic measurements” [l0, and being 
based on a “naive assumption” 03 , concerns no doubt 
based on earlier more explicit criticisms voiced in, for 
example. Refs. 05l. [^. In this paper, we aim to coun¬ 
teract this dismissive attitude to Oklo-based limits on 
Aa = aokio — Qfnow (Henceforth, we shall use a sub¬ 
script 0 to denote a current value, e.g., ao is Ofnow) 

Oklo data constrains shifts AEi = EiOkio — 7?inow in 
neutron capture resonance energies Ei over the interval of 
time since the Oklo natural fission reactors were active 
(about 1.8 Gyr ago). Most attention has been directed 
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at the n + capture resonance nearest threshold. 

Much has been made of the uncertainty associated with 
the modeling of the operation of the Oklo reactors and 
its impact on the values of AE^ extracted; in view of 
the resonance structure of the pertinent neutron absorp¬ 
tion cross-sections, shifts AEi cannot exceed 50meV in 
magnitude 21 |, whereas, on the basis of a representative 
set of reactor model studies (see Table El), we adopt, 
in effect, a 95% confidence level (C.L.) upper bound on 
\AEi\ of 23meV for the ^®°Sm resonance: thus, other 
treatments of reactor dynamics cannot weaken our Oklo- 
inspired bound by more than a factor of 2. We shall also 
revisit the lingering issue of the “second solution” for AEi 
encountered in some analyses of Oklo data and attempt 
to lend further weight to the claim of Ref. that this 
second solution should be ignored. 

As for the reduction of a bound on AEi to a bound on 
Aa, any shift AEi is most appropriately interpreted 
in terms of changes in both a and Xq = mq/A, where mq 
is the average of the u and d current quark masses, and A 
is the mass scale of quantum chromodynamics. Formally, 


and 60, respectively. The values of R and S (with T = i) 
fitted [29j to astrophysical measurements in the direction 
of the radio source PKS1413+135 and atomic clock data 
also suggest that ~ 10 , as do the best fits of the more 
extensive study in Ref. |l7| . but no firm conclusions can 
be drawn until the errors in R and S are substantially re¬ 
duced. The situation is more clear-cut for the many uni¬ 
fication scenarios of Refs, the smallest non-zero 

value of \Rq\ (scenario 6 , 7 = 70) comfortably exceeds 
0.5 (by a factor of more than 3). 

In what follows, we shall assume that Eq. o can be 
used. A lower bound on \ka \ then suffices to establish an 
upper bound on |Aa|/ao- Our focus will be on correc¬ 
tions to the Damour-Dyson estimate of a lower bound to 
I fco, I. We hope to convince even the skeptical reader that 
the Damour-Dyson formula is accurate to better than 
25%. 

The Damour-Dyson lower bound to \k(y\ (in our no¬ 
tation, follows from ingenuous approximations 

which imply the inequality (cf. Eq. (43) of Ref. [l^, in 
which Qi is denoted by Ri) 


AEi = kq 


AXq 

XqO 


Aa 

+ ka - 

Oo 


( 1 ) 


where, for the Sm resonance, it is customary to set kn 
— 1.1 MeV (based on the work of Damour and Dyson [23|, 
who actually conclude that ka < —1.1 ± 0.1 MeV), but, 
unfortunately, kq is poorly known: on the basis of the 
existing estimates ( 2 J,[ 23 , it can be argued that \kq\ ^ 
10 MeV. (It is also conjectured a that kq is approxi¬ 
mately independent of the choice of target nucleus.) For 
the purposes of setting an upper bound on Aa, one can 
work with the related inequality 


|AE,| > \k\R^\ 


l| \ko 


|Aa| 

ao 


( 2 ) 


where k = \kq/kc,\ and i?“ = (AXq/Xqo) / (Aa/ao). 
The more manageable result 


|AE,| > \kc\ 


|Aa| 

ao 


(3) 


will apply whenever = 0 or |i?“| > 2k~^ [so that 
the factor |K|i?q j — 1| on the righthaiid side of Eq. ([5]) 
is greater than or equal to unity]. Given the above re¬ 
striction on \kq\ and our estimate of ka [see Eq. (fTT]) [. 
which together imply k > 4, Eq. Q holds for ^®°Sm if 

\R^\>h 

This criterion can be tested in any realistic model of 
variations in a, which, through its dynamics, must fix 
a particular pattern of correlations between a and other 
fundamental parameters like mq. In the almost u biq - 
uitous phenomenological unification scheme of Ref. [ 23 , 
which can be characterized by the free parameters R, S 
and T, Rq = —R+^{S+1)T. For the canonical choices of 
(i?, 5, T) = (30,160, i) (23 and (36, 240, i) [H], i?“ = 30 


ka < = -■( 4 ) 

where Qi is the equivalent rms radius of the charge distri¬ 
bution of the compound nucleus (CN) state |i) formed by 
thermal neutron capture, and 6g.^{r'^) (> 0) is the differ¬ 
ence between the mean-square charge radii of the ground 
states of the daughter nucleus and of the target nucleus 
(in the case of most interest, ^®°Sm and respec¬ 

tively). As k^^ < 0, Ifcaj is bounded from below by 
if the inequality in Eq. o does, indeed, hold. 

A simple matter, albeit with numerically significant 
consequences, concerns the choice in Ref. j23j of Qi = 
8.11 fm for the ^^°Sm compound nucleus; this value is 
the result of a computation with a formula (Eq. (50) in 
Ref. ( 33 ), which contains not one but two critical tran¬ 
scription errors (from Eqs. (5) and (6) in Ref. [13): a 
more reasonable value, calculated directly from the mea¬ 
sured rms charge radius i?ch [13 of the ground state, 
would be Qi = 6.50 ± 0.01 fm, where the error arises pre¬ 
dominantly from neglect of (5*(r^), defined in the next 
paragraph. [We estimate this contribution to the error 
in Qi with the expression Si,{r‘^) = — («:*o + k*£) implicit 
in Eq. ([9]) , the averages of our results for k*o ) in Ta¬ 
ble Hill and the maximal choice e = 0.05.] From the Sm 
data in Table X of Ref. [13 (which supersedes that em¬ 
ployed in Ref. dl), the difference in mean-square radii 
(5ga(r^) for the isotope pair jg 0.250(20) fm^, 

about 20% larger than the value adopted in Ref. [23| . 
Together, these different parameter values imply that the 
value of k^^ for Sm data should be revised to k^^ = 
-2.51 ±0.20 MeV. 

Our basis for gauging corrections to is an earlier 
inequality [Eq. (36)[ in the analysis of Ref. [13, which 
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relies only on the justified neglect of exchange contribu¬ 
tions to Coulomb energies for its validity. It reads 

ka < J ViSpd^r, (5) 

where Vi is the electrostatic potential of the excited CN 
state |i) and Sp = is the difference between 

the charge densities of |z) and the ground state \t) of the 
target nucleus. In terms of k^^, the right-hand side of 
Eq. (O is 

+ J V,5pd^r + j {Vi-Vu)5pd^r, 

r>Qi 


where Vu is the potential of a uniformly charged sphere 
of radius Qi (and charge Ze), 


V = 


Ze Qi 
Qi r 



and 5*(r^) is the difference between the mean-square 
charge radii of the excited state |i) and the ground state 
of the daughter nucleus. We identify the second through 
fourth terms above as the excitation, Coulomb, and de¬ 
formation corrections to respectively. 

As explained in Ref. [22|, the Coulomb correction 
(i.e., the integral involving Vi) compensates for the use 
in Ref. of the electrostatic potential appropriate to 
the inside of a uniformly charged sphere (of radius Qi) to 
describe the nuclear Coulomb field throughout all space. 
In view of the greater spatial extent of the charge distri¬ 
bution of the compound nuclear state |i), (5*(r^) > 0 and 
the density difference 6p is positive for r > Qi. As a re¬ 
sult, the excitation correction (proportional to —(5*(r^)) 
is negative while the Coulomb correction is positive. The 
deformation correction is found to be positive (see Table 
mi). The consequent partial cancellation of these correc¬ 
tions proves crucial. 

To compute the corrections to k^^, we need densities 
Pk in the vicinity of the nuclear surface {k = n,p,c for 
neutron, proton, and charge, respectively). There is not 
enough experimental data to permit a model independent 
descri ptio n of these densities. Following the example of 
Refs. [SSHSTI, we adopt the deformed Fermi (DF) func¬ 
tions 


Pk = Pok 


1 -I- exp 


r - Ck[l + p 2 kY 2 o{n)] \ 


Zk 


( 6 ) 


Despite the ad hoc empirical origins of this density pro- 
hle, it serves, with some modifications, as a template 
for the extraction from nuclear energy density func¬ 
tional theory of information on surface diffuseness in de¬ 
formed nuclei. Evaluation of the corrections to in¬ 
volves only charge densities, but our scheme for estimat¬ 
ing the charge density parameters of the excited state |i) 


TABLE I. Proton and charge density parameters for 
ground states. 


Isotope 

POp 

(fm-3) 

Cp 

(fm) 

Zp 

(fm) 

P2p 

Cc 

(fm) 

2c 

(fm) 

SkM* 

0.0679 

5.86 

0.502 

0.190 

5.83 

0.560 

SLy4 

0.0679 

5.86 

0.503 

0.184 

5.83 

0.561 

SkM* 

0.0673 

5.88 

0.499 

0.232 

5.84 

0.556 

SLy4 

0.0671 

5.88 

0.501 

0.215 

5.85 

0.559 


presupposes knowledge of neutron and proton densities 
separately. 

Inspection of the DF density parameters obtained in 
experimental studies of even-even Sm isotopes re¬ 

veals, where comparisons are possible (^^°Sm and ^®^Sm), 
some systematic inconsistencies, which suggest that un¬ 
certainties in the surface diffuseness Zc ^-nd quadrupole 
deformation /32c exceed 10% (and could be as much as 
30 % or so). In the analysis of the experiment 

the value of /32c is interpolated from effective /32c values 
for and ^®°Sm (relevant to Coulomb excitation). 

Given these difficulties with empirical DF density param¬ 
eters, we prefer to work with the theoretical Zp’s and /32p’s 
inferred in Ref. from HF-fBCS calculations with a 
contact surface pai ring interaction (constrained by the 
findings of Ref. [39j) and the Skyrme functionals SkM* 
and SLy4. 

Consistent with the restriction to quadrupole defor¬ 
mation in Eq. we disregard the other significa ntly 
smaller deformation parameters determined in Ref. |38l |. 
Their effect on the Coulomb correction to is negli¬ 
gible. The “surface polarization” and the angular depen¬ 
dence of the radial diffuseness identihed in Ref. [38| is 
more of a concern for the other 2 corrections, but these 
features of densities are suppressed by the angular aver¬ 
aging implicit in the calculation of volume integrals. 

As Ref. [s^l deals only with even-even nuclei, we set 
the Zp and /32p parameters for ^^®Sm equal to the aver¬ 
ages of the results for and ^®°Sm (the interpola¬ 

tion scheme of Ref. 37|). We fix the values of Cp and 
pop for both ^^®Sm and by requiring that the pro¬ 

ton density be normalized (to the number of protons) and 
that its second moment {r'^)p reproduce the experimental 
mean-square charge radius, calculated with the standard 
relation (r^)c = {r‘^)p + where the proton rms 

charge radius rp = 0.8775(51) fm and the neutron 
mean-square charge radius = —0.1161(22) fm^ 0. In 
fact, our values of Cp and Zp for ^^°Sm (see Table|l| differ 
only very slightly from those of Ref. [38|, a reflection of 
the size of the deformations we have omitted. 

The charge density parameters Cc and Zc in Table U are 
found under the reasonable assumptions that /32c = P 2 p 
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and pqc = pop] the method of Appendix A in Ref. 42l |. 
generalized to accommodate a non-zero quadrupole pa¬ 
rameter, is used. The SkM* charge parameters for ^®°Sm 
are very similar to the empirical parameters of Ref. . 
(The SLy4 quadrupole deformations for agree 

to within 5% with those of the hnite range droplet 
model ( 4^.1 

The complexity of CN states means that theoretical 
approaches can and must use statistical methods 44l |. 
Recognition of the universal character of properties of 
quantum chaotic systems (of which the compound nu¬ 
cleus is a prototype El) broadens the scope of the ar¬ 
guments that can be brought bear to include insights 
deduced from studies of other more tractable many-body 
systems. Central to our estimates of DF density param¬ 
eters for the state |i) is the widely accepted [46| eigen¬ 
state thermalization hypothesis (ETH) [T^ , according to 
which the expectation value of a few-body observable O 
in an individual state coincides with a microcanonical en¬ 
semble average of O, provided the system manifests fully 
developed many-body quantum chaos. More precisely, 
hnite-size scaling studies [31 indicate that the relative 
error incurred in invoking the ETH decreases as 2 , 
where d is the dimension of the Hilbert space. As the di¬ 
mensionalities d of the many-particle shell model spaces 
needed to describe even low-lying states in ^®°Sm (far 
below neutron threshold) exceed 10® by several orders of 
magnitude, our use of the ETH is justified to accuracies 
of at least 1 part in 10^. Detailed nuclear shell model 
studies (3| hi far smaller model spaces also conhrm the 
appropriateness of the ETH for quantum chaotic nuclear 
states. 

The ETH allows us to adapt the microcanonical en¬ 
semble treatment [IB-liii of mononuclear configurations 
formed in heavy ion reactions to our problem of the densi¬ 
ties pk of |i). The microcanonical analysis can be limited 
to the determination of the surface diffusenesses Zk of |i) 
or the more convenient susceptibilities Xk = Zk/zkg — 1, 
where Zkg denotes the value of Zk for the ^^®Sm ground 
state. Even at excitation energies much higher than 
that of |i), central densities pofe are unchanged [Ulli, 
which means that the equivalent sharp radii Rk, defined 
so that ^RkPko is equal to the volume integral of pk, 
are also unchanged. This assumption, coupled with the 
fact that the quadrupole shape parameters /32fe of If) 
can be constrained by reference to existing studies [541 - 
[3 | of ^^®Sm (see following paragraph), and the relation 
^k = C'Kl + ^/^ 2 fe implies that the central 

radii Ck can be found once the Zk’s are known. 

The effect of excitation on the shape of ^®®Sm (and 
other nuclei) has been studied in thermal relativistic 
mean-field theory 54 - 1^ with the versatile NL3 interac¬ 
tion, large model spaces (with no inert core) and a sound 
description of Coulomb interactions (all improvements on 
earlier investigations). For temperatures up to 0.75 MeV, 


corresponding to an average excitation energy in ^^®Sm 
exceeding the energy of |i), the reported quadrupole de¬ 
formation parameter /32 increases very slightly, due to 
the weakening of pairing correlations, as the tempera¬ 
ture increases (3) S] (/32 denotes the common value 
of /32p and P 2 n)- Similar behavior is observed for ^^"‘Er 
in a finite temperature Hartree-Fock Bogoliubov calcula¬ 
tion js^l, which also uses a realistic effective interaction 
(the DIS Gogny force) and a large conhguration space. 
(The example of ^®^Er furthermore shows that, at these 
temperatures, it is unnecessary to distinguish between 
the mean-held value of (32 and its thermal average: in 
addition to their being numerically close, the trend in the 
mean-held value survives in the thermal average (3|.) 

We interpret these hndings about deformation to mean 
that quadrupole deformation parameters of |f) do not ex¬ 
ceed their values in the ground state of ^®®Sm W more 
than 5%, an estimate based on Fig. 1 in Ref. [^. We 
also take the smallness of thermal huctuations in /32 as 
justihcation for ignoring huctuations in a microcanoni¬ 
cal treatment, i.e., we set the /32fc’s for |i) equal to their 
microcanonical ensemble averages (3 2 k- A® the Coulomb 
correction to displays signihcantly more sensitivity 
to surface diffuseness than to deformation, a more careful 
treatment of huctuations in the Xfe’s is needed. 

Structural features of ^^®Sm inhuencing the values 
of the susceptibilities Xk include the presence of un- 
hlled high-n-low-Z states [3| in the vicinity of the Fermi 
level and surface “tidal wave” excitations [3| compris¬ 
ing rotation-aligned octupole phonons. Excitation of 
the 2.615 MeV octupole vibration in ^®®Pb induces a 
change [giI of only about 1.6% in the surface diffuseness 
of the relevant nuclear potential. To the extent that the 
octupole state in doubly magic ®®®Pb is typical and state- 
dependent changes in a nuclear potential rehect changes 
in the matter distribution, we should not expect the Xfe’s 
for |i) to be more than a few percent or so. 

The energetics of small changes in surface diffuseness 
have been determined (3 | within the self-consistent nu¬ 
clear Thomas-Fermi model. Effects of deformation are 
ignored. Eor \xk\ 0.2, it is found that the change in 
the energy of a nucleus (on Green’s valley of stability) is 
adequately approximated by the quadratic form 

AA = i(18.63MeV)Ai - 2</.2XnXp + , 

where the coefficients (j)i are presented in Table I of 
Ref. [3| as cubic hts in x = (Z/100)^ to numerical re¬ 
sults. 

We require that the most probable values of the Xfc’® at 
any GN excitation energy E* maximize the entropy S. As 
regards the choice of the level density (the logarithm of 
which yields S'), it has been suggested [63| that a compo¬ 
site Gilbert-Gameron (CGG) formula, with a larger en¬ 
ergy shift and a multiplicative enhancement in the Fermi 
gas regime, encapsulates qualitatively the influence of 
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TABLE II. Parameters of a in Eq. Q (flofe’s in MeV and of 
the Gaussian Xp-distributions in the BSFG and GT models. 




Q-On 

K, 



a 

SkM* 

5.44 

6.29 

5.72 

BSFG 

0.00806 

0.0541 

SLy4 

5.44 

6.28 

5.68 

CT 

0.00696 

0.0494 


shell effects, collective excitations and pairing correla¬ 
tions on the level density. In lieu of specific information 
on the magnitude of these modifications for ^^°Sm, we 
employ both the back-shifted Fermi gas (BSFG) and the 
constant temperature (CT) formulas, with param¬ 

eters (appropriate to E* < 10 MeV) taken from Table II 
in Ref. j^. Together, the BSFG and GT models should 
bracket the range of behavior manifested by a modified 
CGC formula. 

For each of these models, we follow Ref. and sub¬ 
tract from E* the diffuseness expansion energy AE: e.g., 
we set the BSFG entropy S'bsfg = ‘2\/a{E* — AE — Ei), 
where Ei is the BSFG energy shift of Ref. (o^l- We also 
incorporate the impact of increasing surface diffuseness 
on the level density parameter a in the BSFG model. 
Guided by the structure of the standard leptodermous 
expansion of a = a„ -I- Op , we demand that 


— ^ aok 
k 


1 + + Afc) 

rifc 


(7) 


where the strength k of the surface terms is fixed so that 
the systematic A-dependence |66| of obsfg is reproduced, 
i.e., Kj2k<^ok^ /J2k^ok = iP/a)A-^ with a and (3 
taken from Eq. (61) in Ref. [^. The oofe’s are chosen so 
that ao^op = {N/Z){pQp/pQn)i (cf. Eqs. (13) and (14) 
in R ef. 1651 1. and Eq. ([7]) reduces to the value of obsfg in 
Ref. [6JI when the Xk’s are zero. [Instead of a single fac¬ 
tor K, Eq. © should, in principle, contain two deforma¬ 
tion dependent factors Kk, approximately proportional to 
1 + but the ratio Kpjun differs from unity by less 

than 0.5%.] Einally, we employ the relation between 
Obsfg and the temperature parameter in the CT model 
to substitute this temperature parameter in the denomi¬ 
nator of ^CT by 5.164a“°'^®^, where a is given by Eq. 0. 
As the lefthand half of Table m illustrates, our values of 
the flofe’s and k are insensitive to whether we adopt the 
SkM* or the SLy4 ground state densities of Ref. [S^ ■ 

Eor each choice of S, the related microcanonical prob¬ 
ability distribution function for the Xfe’s is well repre¬ 
sented about its maximum by a bivariate Gaussian. The 
means p and standard deviations cr of the associated 
Gaussian marginal distributions for Xp, listed in Table HIl 
are the same (to 3 significant figures) for our two sets of 

(uop, UOn, ^)• 

We can now discuss estimates of the corrections to 
We begin with the excitation correction, propor¬ 
tional to Si,{r^). The relation between (r^) = (r^)c and 


{r^)p quoted above allows us to set 5i,{r^) = 5i,{r‘^)p. 
There are two contributions to 6i,{r'^)p, one proportional 
to (5*/3|p, the other proportional to i5*Zp. In line with our 
earlier assumptions about /32p, we approximate (5*/3|p as 
Pip ~ P 2 pg^ where (32pg is the value of /32p for the ^®°Sm 
ground state. We calculate <5*Zp by averaging over the 
Gaussian distribution for Xp- Thus, 


= {-Rl 


y 


2p 



3-T^'^Zpg{‘^P +pP + ( 8 ) 


For the small range of /32p values we admit (from /32pp to 
1 .05,52pg)j the excitation correction is excellently approx¬ 
imated by a linear function of e = P 2 p/P'Zpg ~ which it 
is natural to write as 


(Ze)^ / / X 

++ K*^)- (9) 

The constant coefficients k*o Eind are given in Ta¬ 
ble HIl 

Our calculations reveal that the Coulomb and deforma¬ 
tion corrections can also be regarded as linear functions 
of e. For ease of comparison, we adopt parametrizations 
of the same form as in Eq. (j9]) with k*-o , replaced by co¬ 
efficients Kio, k' {i = c,d), again tabulated in Table Hill 
[Here, c (d) denotes the Coulomb (deformation) correc¬ 
tion.] 

Concerning the Coulomb correction to the inte¬ 

gral I of the product of the function Vi with a DE charge 
density over the volume outside a sphere of radius 
Qi can be reduced to the angular average of a linear com¬ 
bination of complete Fermi-Dirac integrals Tkiv) [13 ■ 


jzer 

3Q^ 



where pa = [Qi - Cca[l + P 2 caY 2 o[Vt)\) /Zca and ■ )n 
denotes the average over all solid angles H. The Coulomb 
correction is 5^^! -\- 6i,I. The dg^J-term is computed with 
the charge distribution parameters of Table HI In the d*/- 
term, which entails averaging over the Xc-distribution, 
we substitute /32ci by the range of /32p-values considered 
above and assume that the Xc-distribution can be approx¬ 
imated by that of Xp- The Xc-dependence of the central 
radius Cd is also taken into account. 

Calculation of the deformation correction is facilitated 
by multipole expansions in even I spherical harmonics 
Yiq. For convergence to 3 significant figures, the Z = 0, 2, 
and 4 terms suffice. The quadrupole contribution to the 
deformation correction is dominant. The Xc and f32ci 
dependence is dealt with in the same way as that of the 
Coulomb correction. 

Together, our results for ^®°Sm in Table ITTTI imply that, 
depending on the actual value of e and the choice of 
model (SkM*+BSFG, etc.), the excitation and deforma¬ 
tion corrections are separately somewhere between 15% 
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TABLE III. Individual correction coefficients ACio, (* = c, *, d) for data and net correction coefficients kq, k' for ^^*^ 8111 , 

^®®Gd and ^®®Gd data. (All coefficients are in units of fm^.) 






150 

8m 





156 

Gd 

158 

Gd 



^cO 

n'c 

«;*o 

f 

KdO 

/ 

Ko 

k' 


k' 

/^o 

K 

8kM* 

B8FG 

0.0320 

0.119 

-0.0469 

-0.964 

0.0672 

0.543 

0.052 

-0.30 

-0.005 

- 0.68 

-0.018 

-0.75 


CT 

0.0302 


-0.0403 


0.0669 


0.057 


0.000 

M 

-0.014 

M 

8Ly4 

B8FG 

0.0305 

0.101 

-0.0474 

-0.833 

0.0400 

0.452 

0.023 

-0.28 

0.012 

-0.63 

-0.016 

-0.72 


CT 

0.0287 


-0.0407 


0.0397 


0.028 


0.016 

M 

-0.011 

M 


and 40% of the revised value of and the Coulomb 
correction is about lO-to-15%. However, because of the 
partial cancellation of these corrections, the net correc¬ 
tion A:™''*' is relatively modest: fc™"'’' is always less than 
25% of the revised value of Following the inter¬ 

model strategy of Ref. [g^ for the estimation of the error 
associated with the use of models, we compute the net 
correction for the two extreme values of e (i.e., 0 and 
0.05) in all four of the models in Table ITTTl and employ the 
average of these values and their standard deviation as, 
respectively, our best estimate for and the related 
uncertainty: 


= 0.33±0.16MeV. 


( 10 ) 


k„ < k^'^ = fcDD + 


,DD 


(More robust statistics |69|, namely, the median and the 
median average deviation from the median, yield almost 
the same numerical result: fc™’'"' = 0.33 ± 0.14MeV.) 

The preceding analysis implies that, for the ^^°Sni res¬ 
onance of interest, 

=-2.18 ±0.26 MeV, (11) 

where the errors in k^^ and A:™'*’ have been added in 
quadrature. 

The generic character of our arguments about correc¬ 
tions to k^^ means that they can be adapted to other 
complex nuclei, in particular, the well-deformed ^^®Gd 
and ^®®Gd isotopes considered in Ref. [^. Using the Gd 
data in Refs, (sj, 34|, k^^ = —1.1 ± 0.1 MeV (—1.3 ± 
0.2 MeV) for thermal neutron capture by ^^®Gd (^®^Gd). 
The differences in A:)^^-values are primarily a conse¬ 
quence of the variable extent of odd-even staggering in 
mean-square radii, an effect which cannot be reproduced 
by an approximation based on the Coulomb term in 
the Bethe-Weizsacker mass formula. The only aspect of 
our treatment of corrections which requires modification 
is the choice of maximal value for e. As the ground state 
quadrupole deformations P 2 of ^®®Gd and ^®®Gd are very 
close to that of we appeal to the thermal behavior 

of (32 for in Fig. 2(a) of Ref. to constrain e 

to the interval between 0 and 0.03. The net corrections 
to the fc^^-values above are then, parallelling the ear¬ 
lier Sm analysis, —0.04±0.13MeV and —0.26±0.11 MeV 
for ^®®Gd and ^®®Gd, respectively. The corresponding 


limits on are ka < —1.1 ± 0.2 MeV for ^®®Gd and 
ka < -1.6 ± 0.2 MeV for issQd. 

These results for Gd isotopes have a bearing on the 
resolution of the “second solution” problem Q . From the 
comparison (7ll | of measured rare earth isotope abun¬ 
dances in Oklo sample KN50-3548 with calculations 
based on present-day neutron absorption cross sections, 
it can be inferred [2l| that shifts AEi in resonances must 
be less than 50 meV in magnitude. Nevertheless, the re¬ 
sults of Damour and ERson for Oklo samarium data 
alone can be interpreted [8| to mean that AEi lies in 
either a “right branch” AEi = 46 ± 22 meV, compatible 
with zero, or a “left branch” AEi = —94 ± 13 meV, incon¬ 
sistent with the 50meV bound on |Ai5i|. Subsequently, 
for Oklo RZIO and RZ13 samples, a similar non-null so¬ 
lution was ^ain found (in addition to a null solution) W 
Fujii et al. [3| [AEi = —97 ± 8 meV) and Gould et al. Q 
{AE, = -90.8 ± 5.6 meV). 

Exciting as the prospect of a non-zero result for AEi 
may be, we are of the opinion that the intervals which 
do not overlap with zero are an artifact, ultimately, of 
the symmetry of a Breit-Wigner absorption cross section 
about the resonance energy. If there is a (physical) energy 
interval to one side of a resonance which can be associated 
with a particular range of effective capture cross section 
values, then there will necessarily also be an unphysical 
interval on the other side of the resonance. 

An attempt to reconcile Oklo data on thermal neutron 
capture by all three of the isotopes ^^^Gd, and 

^®^Gd reinforces this point. Generalizing the analysis of 
Ref. 0 to include the AVg-term in Eq. o, one expects 
that AE| = AEi — kaAa/ao should be approximately 


TABLE IV. Icr intervals for AEi and their sources. 



Ei 

Aa/ao 

AEi 



(meV) 

(meV) 

(meV) 

Ref. 




4± 16 

18J 

n + i«8m 

97.3 

2.2 

7.2 ±9.4 

[9] 




1.9 ±4.5 


n-f i®®Gd 

26.8 

0.9 

-8.5 ±17.5 

[8] 

n + i®"Gd 

31.4 

1.1 

-8.5 ±17.5 

18J 
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the same for all nuclei 2^. For the sake of argument, 
we calculate the Aa-contribution to using a quasar- 
based estimate of Aa ~ —8.6 x 10“^°anow uM the 
values above of (In the quasar-based estimate, Aa 
is attributed to the motion of our local galaxy cluster 
relative to the Australian dipole in the time since the 
Oklo reactors were active.) The corresponding values 
of the Aa-term, presented in Table IlYl are an order of 
magnitude smaller than the AEfs used in Ref. Q (also 
given in Table El). Hence, the conclusions in Ref. @ 
about AEfs will apply to the AE^’s: if one admits the 
presence of post-reactor contamination in the Gd data 
(at the 3-to-4% level), then one can isolate A£’)'-intervals 
for all three nuclei which are approximately the same in 
as much as they all overlap zero, whereas the unphysical 
interval is negative for the case of Sm and positive for the 
Gd isotopes. This pattern will continue to apply if the 
actual values of the Aa-term are used provided, of course, 
that they do not differ substantially from the choices in 

Table E 

Despite the many uncertainties to which the analysis 
of Oklo data is subject, the different Sm results for AEi 
in TableJIVI agree to within a factor of 2 with the result 
of Ref. 1^, which we adopt. If we combine the bound 
of Ref. [9fl on AEi with our restriction in Eq. HH) on 
ka, and simplify the distribution of the ratio — AEj. / k^‘^ 
along the lines of Sec. 4 of Ref. [^, we deduce the upper 
bound at 95% G.L. of 


|Ac 


CkQ 


< 1.1 X 10” 


( 12 ) 


which is comparable to the Oklo-based limits listed in 
Refs. and [t^, but on a sounder footing. The quasar- 
based prediction in Ref. [73 of |Aa|/ao is compatible 


with Eq. m- 

Assuming a linear time dependence for a over the last 
1.8 billion years, Eq. m implies that the present-day 
time variation of a is subject to the (95% C.L.) bound 


da 

dt 


< 0.61 X 10"^^yr-\ 

0 


which is an improvement on the best limit from atomic 
clock experiments [l^ . 

In this paper, we have been at pains to demonstrate 
that the order of magnitude of the bound in Eq. m 
is reliable. We believe that neglect of the Oklo-based 
bound on Aa is unfortunate. It provides a restrictive 
low-z datum which can help to select from the current 
plurality of theoretical models admitting variations in a, 
those which are phenomenologically acceptable. Most of 
the model studies which have included the Oklo limit on 
Aa in their analysis, have been content to invoke the 
result of |Aa|/ao 10“^ to be found in Damour and 
Dyson’s seminal paper [i^. It would be interesting to 
see how previous conclusions are revised if a bound on 


|Aa|/an of the order of 10 ® is adopted. Models [hi . 


75l . l76| which naturally suppress the variation of a in the 


presence of matter may well be preferred to the exclusion 
of all others. It should also be instructive to consider the 
impact of this bound on feasibility studies pertaining to 
the ambitious program of astrophysical measurements of 
the redshift dependence of parameters like a put forward 
in Ref. 771 and reviewed recently in Ref. Q. 
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